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A detailed study of the lattice gas (Ising) model shows that the correct definition of cluster 
concentration in the vapor at coexistence with the liquid leads to Fisher plots that are unique, 
independent of fixed particle number density pfi xe d, and completely reliable for the characterization 
of the liquid- vapor phase diagram of any van der Waals systems including nuclear matter. 
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The liquid to vapor phase transition in nuclei and in 
nuclear matter is a time honored topic nearing its cul- 
mination with a complete phase diagram of bulk nu- 
clear matter Q. On one hand, general indicators of 
phase transitions in finite systems have been explored 
theoretically @, H H @, while, on the other, experi- 
mental evidence for phase transitions has been claimed 
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Specifically, one method based on Fisher's droplet 
model | l2| was applied to experimental nuclear multi- 
fragmentation data and the full liquid-vapor coexistence 
line (phase diagram) from low temperatures up to the 
critical point was extracted Q, |n| . 

The validity of such an approach has been demon- 
strated by applying it to the cluster concentrations gen- 
erated by the Ising model The cluster concentra- 
tions collapse onto a sing le line (a Fisher plot) accord- 
ing to Fisher's theory and allow for the quantitative 
recovery of the Ising phase diagram and the associated 
parameters, e.g. the critical temperature T c . 

This success has been tempered by the claim that 
the Ising cluster concentrations, while giving excellent 
Fisher plots, depend upon the mean particle number 
densities pfi xe d for which they have been calculated, and 
project different critical points for different mean densi- 
ties |3 • The implication of these claims is that the 
phase diagram for the Ising model and for nuclear mat- 
ter obtained by such procedures is not unique and hence 
incorrect. 

In this Rapid Communication it will be shown that 
such claims arise from an improper use of Fisher's theory, 
and that a correct application of that theory leads to a 
unique Fisher plot, independent of the value of pfi xe d, 
which therefore can be directly transcribed into a unique 
phase diagram with the correct critical parameters. 

Preliminary to the present analysis of cluster concen- 
trations of the Ising model, it is worth recalling that, for 
an infinite system, thermodynamics demands that the 
saturated vapor properties be independent of the mean 
density of the system. Consider a cylinder with a given 
amount of water and a piston that can alter the volume of 
the fluid. At any fixed temperature T below the critical 
temperature, the vapor in equilibrium with the liquid has 
a pressure p(T) which is completely determined by the 
temperature (univariance) . The position of the piston, to 
which the appropriate counter pressure p(T) is applied, 



is in indifferent equilibrium (and the system has an infi- 
nite isothermal compressibility), i.e. it can be moved up 
or down without experiencing any restoring force, thus 
covering the entire density range from liquid pi to vapor 
p v until the last drop of liquid remains. The properties of 
the vapor in the univariant regime depend exclusively on 
T and not on the total amount of vapor or liquid in the 
system because at phase equilibrium (coexistence) the 
chemical potential is uniquely defined by T. Thus the 
corresponding Fisher plot depends only on the property 
of the vapor and must be independent of Pfixed- 

After the last drop of liquid disappears the vapor is free 
from such a constraint. As the piston further rises after 
the disappearance of the liquid, the system becomes bi- 
variant, and the piston is no longer at indifferent equilib- 
rium, but must obey the equation of state: p = /(p, T). 
For the ideal gas limit p oc pT . In these conditions, de- 
viations from the Fisher plot are to be expected unless 
the bivariance is accounted for in terms of the chemical 
potential. Thus the obvious caveat: the standard Fisher 
plot is unique and independent of the mean density pro- 
vided only cluster concentrations at coexistence are con- 
sidered. 

The features discussed above can be observed in clus- 
ter concentrations generated from fixed magnetization 
(Mfi xe( i) Ising calculations on a two dimensional (d — 2) 
square lattice of side L. The calculations are performed 
via the method of Heringa and Blote [pfij. Periodic 
boundary conditions and a lattice side of L = 80 are 
used to minimize the finite size effects of the container; 
a shift of only ~ 0.45% in the value of T c is expected 
for this system [T?l li"8| . The ground state liquid drops, 
A d (T = 0) = 1280, 640 and 320 give fixed mean densities 
(Pfixed = A d (T = 0)/L d ) Of Vfixed = 0.2, 2 p fixod = 0.1 
and 3 pfi xe d = 0.05, respectively. These drops are large 
enough that their finite size effects are small, though still 
observable. Finite size effects associated with the drop 
size both in the Ising model and in the general discus- 
sion are neglected here in order not to cloud the main 
issue of this paper. However, they can be dealt with by 
employing the complement approach which is discussed 
elsewhere 0. 

In these calculations the lattice acts as the container 
in which a liquid drop is placed. The fluid is made up 
of the minority spins (e.g. down spins) with the vacuum 
being the majority spins (e.g. up spins). At the lowest 
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FIG. 1: Temperature T as a function of molar volume of 
liquid (left) and vapor (right) for different fixed mean densities 
Pflxcd compared to the analytic result. Deviations from the 
analytic result are due to finite size effects of the drop. 



FIG. 2: Pressure p as a function temperature T. The pres- 
sure stays on the coexistence curve until complete evaporation 
and ideal behavior is observed. Deviations from the analytic 
solution are due to the finite size of the drop. 



temperatures the up spins congregate into a single clus- 
ter, the liquid drop, surrounded by vacuum. At higher 
temperatures the vacuum is filled with a vapor of clus- 
ters of down spins, and the liquid drop diminishes in size 
accordingly. Clusters in the vapor are identified via the 
Coniglio-Klein (CK) algorithm 19] to insure that they 
are physical (i.e. the cluster concentrations return Ising 
rather than percolation critical exponents). The largest 
cluster is identified as the liquid drop. It is identified 
geometrically (all like spin nearest neighbors bonded) in 
order to capture the skin thickness associated with liq- 
uid drops [l|. Alternatively, the CK algorithm can be 
applied to the large cluster as well, but the resulting ex- 
cess concentration of small clusters associated with the 
surface of the drop should not be used in evaluating the 
vapor properties. 

For every temperature T and magnetization Mfl xe( j, 
10 5 thermalized realizations were generated to produce 
the cluster concentrations analyzed here. Temperatures 
are in units of spin-spin interaction strength j, densities 
are in units of spins per lattice site, molar "volumes" 
(more properly molar surfaces since the calculations are 
two dimensional) are given in terms of lattice sites per 
spin, pressure is given in units of j per lattice site and 
cluster sizes A are the number of spins in the cluster. The 
standard translation from Ising model to lattice gas is 
made through the density-magnetization relation pfi xc d = 
(1 - M fixed ) /2. 



The density of the vapor p v is 

A 

where ua(T) is the concentration of clusters of size A in 
the vapor. More specifically ua{T) is the yield of vapor 
clusters normalized to the area available to the vapor (or 
volume in three dimensions) and not the total area of the 
lattice. The area available to the vapor is the area of the 
lattice L 2 minus the area taken up by the liquid drop: 
Ad(T) + A^f ' A&(T) + nhoics, where the second term ac- 
counts for the skin thickness of the liquid drop and the 
third term takes into account holes (bubbles) in the liquid 
drop. The molar volume of the vapor is V£ = l/p v . The 
density of the liquid drop is pi — A d (T)/ (A d (T) + nhoies) 
and the molar volume of the liquid is = 1/pi- The 
exclusion of the surface (4-^/ Aa(T)) of the drop in deter- 
mining its density is in keeping with the idea that is has 
some skin thickness. 

The pressure of the vapor p is determined by assuming, 
as in physical cluster models, that all interactions are 
exhausted by cluster formation so that the clusters can 
be treated as an ideal gas and 

p = T$> A (T). (2) 

A 

Performing d = 2 calculations has the benefit that the 
V^i" and p values obtained from the cluster concentra- 
tions can be compared to the analytical result given by 
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Lee and Yang u sing Onsager's solution for the d = 2 
Ising model [20l l'2ll I22I |23| . Figures Q and |21 present 
such comparisons. 



- Analytic solution 
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FIG. 3: The specific heat CV (top), liquid drop size nor- 
malized to the ground state A d (T)/A d (T = 0) (middle) and 
fluctuations in A d (T)/A d (T = 0) (bottom) against T. Peaks 
in Cv and |A (A d (T)/A d (T = 0)) /AT coincide with the de- 
viation from coexistence. Dotted lines guide the eye. 

Figure ^ shows phase separation at low temperatures 
for the three pfi xe d calculations: there is a vapor (right- 
most open symbols) , with the same molar volume as that 
given by the analytic solution (solid curve), coexisting 
with a liquid drop, which has the molar volume of the liq- 
uid as given by the analytic solution. As the temperature 
increases, the difference between molar volumes of liquid 
and vapor decreases, but, because of the constraint im- 
posed by /Jfixed, at some point leaves the coexistence 
curve and attains a constant value equal to l/pfi xc d- The 
smaller the ground state drop, i.e. the smaller the total 
number of particles in the system, the lower the tem- 
perature at which the system leaves coexistence. Devi- 
ations from the analytical solution of V m (T) before the 
departure from coexistence are due to the finite size of 
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FIG. 4: Vapor concentrations ua(T) at different temperatures 
T for each pfi xe d calculation. The curves show ua(T) (heavy 
solid) and ua{T c ) (light dashed) from Eq. 10. The symbol 
definitions are the same as in Fig. [3] 



the drop, which become pronounced at large temperature 
where the drop is smallest. The deviations increase with 
increasing temperature because the drop size decreases 
with increasing temperature (as shown in Fig. thus 
increasing the drop's finite size effects. These finite size 
effects can be accounted for in terms of the complement 
PJ which reduces to the Gibbs-Thomson correction [24| 
for large drops. 

Figure [21 shows that at low temperatures all three sys- 
tems agree with the analytical T-p coexistence curve. At 
high temperatures all three systems develop nearly ideal 
gas behavior and p oc pt^ ed T as shown by dashed lines. 
Each system leaves the coexistence line when the drop 
has evaporated completely; this happens at a higher tem- 
perature for the larger ground state drops. As usual, de- 
viations from the analytical solution of p(T) arise from 
the finite size of the drop. 

The departure of the pfi xe d calculations from coexis- 
tence gives rise to a peak in the fluctuations of various 
quantities which is not related to criticality. In fact the 
peak in the specific heat in pfi xe d d = 3 lattice gas cal- 
culations has been used to map the coexistence curve in 
the p-T plane [2{|. Indeed, the top panel of Fig.[3]shows 
that the peak in the specific heat CV occurs just as the 
system leaves coexistence. The bottom panel of Fig. 
shows that this coincides with the peak in the fluctua- 
tions of Ad(T). The peaks occur somewhat before the 
Pfi X od constraints are felt because of the effects of the fi- 
nite size of the liquid drop. Evidence for this can be seen 
by comparing the peak locations in Fig. [3] with the depar- 
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FIG. 5: [color online] Fisher plots: the scaled cluster concentrations fit with Eq. © (solid line). 



ture from coexistence in Fig. ^ The coincidence is more 
precise as A^(T = 0) increases. The peak in fluctuations 
observed here do not arise from criticality but because 
coexistence has ended. 

Figure 0] shows the vapor concentrations of the three 
systems and the prediction of Fisher's cluster model 0] 
for vapor concentrations at coexistence: 

n A (T) = q A- T exp(-^^j (3) 

where qo is the normalization, r is a topological criti- 
cal exponent, c is the surface energy coefficient, a is 
the perimeter to area exponent (or surface to volume in 
d = 3) and e = (T c - T)/T c . The solid curves in Fig. H 
are obtained from Eq. J3J using the d = 2 standard pa- 
rameter values given in the first column of Table [I] and 
Co = 8. The dashed curves in Fig. 0] show the critical 
point power law. 

Figure 0] shows the main point of this paper: the va- 
por cluster concentrations are well reproduced by Eq. 
when the system is at coexistence, but not when the sys- 
tem is away from coexistence; in other words, the location 
of the pfixed-T 1 point in the p-T plane affects the vapor 
concentrations. In Fig. 0Ji all three systems are at co- 
existence. All manifest the same vapor concentrations 
described by Eq. ©. In Fig. the 3 pfi xc d system is 
departing coexistence (and entering the large fluctuation 
region) while Vfixed and 2 /?fixed remain at coexistence. 
Accordingly, 1 pfi xo d and 2 /9fi xc d manifest the same vapor 
concentrations described by Eq. @ while 3 /9fi xo d deviates 
above. In Fig.^J:, 3 /9fi xe d is well out of coexistence, 2 /9fi xe d 
is just leaving coexistence. Accordingly, only Vfixed is 
still at coexistence and manifests the vapor concentra- 
tion of Eq. In Fig. QJl, both 2 pfi xe d and 3 /9fi xc d are 
well outside the coexistence region while Vfixed, leaving 
coexistence near T c shows a vapor concentration close to 
the critical point power law. 

In order to further verify that the coexistence vapor 



concentrations are well described by Eq. J3J), the clus- 
ter concentration of each pfi xc d calculation was indepen- 
dently fit to Eq @ leaving cq and T c as fit parame- 
ters, using qo = \/2Q(t — 1) and fixing r and a from: 
(r - 2)1 a = = 1/8 and (3 - r)/a = 7 = 7/4. Clusters 
of size 7 < A < 21 were included in the fit. In the ther- 
modynamic limit the highest temperature that can be 
described by Eq. © is the temperature at which the sys- 
tem leaves coexistence: T « 2.2, 2.05 and 1.6 for 1 pfi xe d, 
2 pRxed and 3 pfl xo d, respectively. However, these systems 
are not at the thermodynamic limit and the fluctuations 
grow large before those temperatures, therefore the high- 
est temperatures fit were T = 1.7, 1.6 and 1.5 for Vfixedj 
2 Pfi XG d and 3 p fixod , respectively. 

TABLE I: Parameter values and fit results 

Pfixod free 02 OA O05 

xl - 1-85 1.41 1.56 

c - 8.7 ±0.1 8.2 ±0.2 7.8 ± 0.2 

T c 2.26915 2.24 ±0.01 2.25 ± 0.02 2.25 ± 0.03 

p c 0.1097 0.114 ±0.003 0.115 ± 0.002 0.12 ±0.01 

C F 0.0966 0.1021 ±0.0007 0.102 ± 0.001 0.102 ±0.001 

The results are shown in Fig.0and summarized in Ta- 
ble[I] All systems show a good collapse to a single Fisher 
plot. All systems return, to within error bars, the same 
T c value close to the the Onsager value, also given in 
Table 13 This is at variance with reference 0]> which re- 
ports changes in T c by as much as ~ 50% with changes in 
mean densities. The present results show, however, that 
by selecting vapor concentrations (given by a proper def- 
inition of vapor and liquid and the correct measure of the 
volume available to the vapor) at coexistence a consistent 
value of T c for all mean densities can be obtained. The 
change in cq observed between the systems is due to the 
finite size of the drop and can be accounted for by the 
complement approach Q]. 
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FIG. 6: The T-p (top) and T-p (bottom) coexistence curves. 
The solid curves show the analytic solution, symbols show the 
results of summing Eq. with the fit parameters in Table[I] 



Figure El also shows the results for clusters of size 
A = 1-6, though they were not included in the fits. These 
clusters are excluded from the fits because their perime- 
ter (or surface in d = 3) energies are poorly described 
by Fisher's ansatz: -E pcr imctcr = cqA° . For example, 
monomers have a perimeter of 4, and (with the d = 2 
Ising perimeter tension of 2) -Ep er i meter = 8; dimers have 
a perimeter of 6, and -Eperimeter = 12; trimers have a 



perimeter of 8, and Eperimeter = 16. Clusters of A = 4 
can have a perimeter of 8 or 10 and in the temperature 
range considered have an average perimeter of A ~ 9, 
and -Eperimeter ~ 18. For clusters 1 < A < 4 these values 
of -Eperimeter were used in place of CoA a and for clusters 
of size A = 5 and 6 -Eperimeter = coA a was used. 

The fit results can then be used to determine the re- 
duced phase diagram accordingly: 
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and 
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(5) 



The sums above and those used to determine qo from r 
were evaluated via methods similar to those in |2q| . 

To determine the absolute density p the reduced den- 
sity from Eq. (0J was multiplied by the known criti- 
cal density of the system p c = 1/2. The top plot in 
Fig. [S] shows the absolute vapor density at coexistence 
as a function of temperature. To obtain the bottom 
plot of Fig. H3 the T-p coexistence curve, the values 
of T c and p c are used with the compressibility factor 
C f = Pc/PcT c = Y,A n A{T)IY,A An A{T)- All systems 
return the same coexistence curve and agree with the 
analytical solution to within error bars. 

In conclusion, it has been shown that: in agreement 
with thermodynamic expectations, the properties of the 
d = 2 lattice gas (Ising) vapor at coexistence are indepen- 
dent of mean particle number density; their Fisher plots 
are unique; and lead to unique critical temperatures and 
phase diagrams. 
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